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We consider a mixture of two species of spin-1 atoms with both interspecies and intraspecies spin 
exchanges in a weak magnetic field. Under the usual single mode approximation, it can be reduced 
to a model of coupled giant spins. We find most of its ground states. This is a complicated problem 
of energy minimization, with three quantum variables under constraints, i.e. the total spin of each 
species and the total spin of the whole mixture, as well as four parameters, including intraspecies 
and interspecies spin coupling strengths and the magnetic field. The quantum phase diagram is 
^ | very rich. Compared with the case without a magnetic field, the ground states are modified by a 

magnetic field, which also modifies the ground state boundaries or introduces new crossover regimes 
*y • on the phase diagram. Without interspecies spin coupling, the quantum phase transitions existing 

in absence of a magnetic field disappear when a magnetic field is applied, which leads to crossover 
regimes in the phase diagram. Under ferromagnetic interspecies spin coupling, the ground states 
remain disentangled no matter whether there is a magnetic field. For antiferromagnetic interspecies 
spin coupling, a magnetic field entangles the ground states in some parameter regimes. When the 
^ ' , intraspecies spin couplings are both ferromagnetic, the quantum phase transition between antifer- 

romagnetic and zero interspecies spin couplings survives the magnetic field. When the intraspecies 
bJT)' spin couplings are both antiferromagnetic, a magnetic field induces new quantum phase transitions 

between antiferromagnetic and zero interspecies spin couplings. 

C$ ' PACS numbers: 03.75.Mn, 03.75.Gg 
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I. INTRODUCTION 



Spinor Bose gases have been extensively studied since a decade ago when it was discovered that they display 
remarkable spin correlations because of spin-exchange scattering between atoms [l|-@, lij-fllj]. However, there have 
^ | not been many investigations on many-body phenomena in mixtures of different spinor Bose gases with interspecies 
spin exchanges. In the first instance, spin-exchange scattering between distinguishable atoms have been less studied, 
perhaps because of incomplete information on inter-atomic potential. To motivate more interests in this direction, 
we note that interspecies spin-exchange interaction can be si gnifi cant. The previous experiments on multi-component 
t^J- ' Bose gases often had atom loss due to spin exchanges [ill Il2j. There were early calculations indicating that the 
, cross-sections of spin-exchange scattering between different atoms may not be smaller than those between identical 
atoms [13j . Recently, there were more calculations as mentioned in the following, though motivated by studying a 
mixture of two species of atoms in frozen spin states, for which spin-exchange scattering is regarded as inelastic. A 
calculation for 23 Na- 85 Rb scattering indicated a quite large difference between scattering lengths of electronic singlet 
and triplet states [3] . This paper also reported a small value of such a singlet-triplet difference of scattering lengths 
for 23 Na- 87 Rb scattering, but contrary result was later reported [l^]. Another calculation found significant singlet- 
triplet differences of scattering lengths for X- 133 Cs scattering, where X= 6 Li, 7 Li, 39 K, 41 K, 85 Rb and 87 Rb |lr|. 
Experimentally, significant differences between singlet and triplet scattering lengths have been observed in 41 K- 
87 Rb, A0 K- 87 Kb and 6 Li( 7 Li)- 23 Na mixtures [l7l - l20j . implying significant interspecies spin exchanges. Spin-changing 
scattering was also observed in 7 Li- 133 Cs [2l| . Moreover, heteronuclear Feshbach resonances can be implemented |22| . 
which can enhance both elastic and inelastic collision rates (23l - [25j . To our understanding, some recent experimental 
set-ups on multi-species Bose gases have come to close to what we need to realize a mixture of two species of spinor 
gases with interspecies spin exchange (26[. 

Further researches on mixtures of spinor gases with interspecies spin exchanges can be motivated by novel many- 
body quantum phenomena in such mixtures, as demonstrated first in a model of a mixture of pseudospin-i atomic 
gases, where interspecies spin exchange leads to richer ground states and phenomena, especially Bosc-Einstein con- 
densation (BEC) with interspecies quantum entanglement, which was dubbed entangled Bosc-Einstein condensation 
(EBEC) |27H30| . As usual, two subsystems constituting the total system are entangled if the state of the total system 
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is not a direct product of those of the subsystems. Otherwise, they are called disentangled. 

This line of researches has been extended to a mixture of two species of spin-1 atomic gases JUHl], in which the 
interspecies spin coupling is simply of Heisenberg form 33] . In the usual approach of single orbit- mode approximation, 
most of the exact ground states in absence of a magnetic field have been found [3lL [33]. However, in presence of a 
magnetic field, only two special parameter regimes have been considered [3l| . Given that the magnetic field effect is an 
important issue, in the present paper, we systematically study the ground states in presence of a weak magnetic field, 
and find out how a magnetic field affects the ground states and phase diagrams of a spinor mixture with interspecies 
spin exchanges. 

The rest of the paper is organized as the following. To make the paper self-contained, we set the stage in Sec. UU 
Then we discuss in Sec. IIIII the ground states of a mixture of two spin-1 Bose gases in a magnetic field, but without 
interspecies spin coupling. In Sec. IIV1 we find the ground states of a mixture with ferromagnetic interspecies spin 
coupling in a magnetic field, based on the calculations detailed in the Appendix A. For antiferromagnetic interspecies 
spin coupling c ab > 0, we divide the range of c ab to three intervals. In Sec.fV] based on the calculations detailed in the 
Appendices B, C, D and E, we find the ground states of a mixture with < c ab < 2jB, where 7 is the gyromagnetic 
ratio, B is the magnitude of the field. In Sec. IVI1 we make some brief discussions on the regime c ab > 2jB. In 
Sec. IVIll quantum phase transitions are described. The issue of characterizing interspecies entanglement is discussed 
in Sec. IVIII1 A summary is made in Sec. IIXI 



II. THE SYSTEM 



Consider a mixture of two species a and b of spin-1 atoms, whose numbers N a and N b are conserved respectively. 
The single-atom Hamiltonian of species a (a — a, b) is 

h 2 

h a = - — \7 2 + U a (r)~ la B-F a , (1) 

where B is a uniform magnetic field, m a , 7 Q > 0, U a and F Q are the mass, the gyromagnetic ratio, the external 
potential and the single-spin operator, respectively, for an atom of species a. With ij) ap representing the field operator 
corresponding to spin \x component of species a (/z = — 1,0,1), the many-body Hamiltonian is 

H= U a + H ab , (2) 

a—a^b 

where 

H a = J dnp^hair)^^ + ^ J drip ap ip ap (c^6 pu 5 pcr + ^F afJ , u ■ F apa )t/j a(T ip av (3) 
is the usual Hamiltonian of spin-1 atoms [l[, 

H ab = j dr^ a ^l p (cfS^8 pa + c ab F ailv ■ F bpa ^ ba ^ av (4) 



is the interspecies interaction (3l| , where Cq and are expansion coefficients in terms of powers of dot product of the 
single-spin matrices of two atoms of species a and are linear combinations of singlet and triplet scattering lengths, 



ng lengths, 



c-2 is proportional to the differences between triplet and single scattering lengths of intraspecies scattering 
and (S^ are similar quantities for scattering between an a-atom and a 6-atom, c?, b is proportional to the differences 
between triplet and single scattering lengths of interspecies scattering, and it has been shown that the coefficient of 
(F n -F b ) 2 is zero 0. 

For each species and each spin state, we follow the usual single mode approximation for the single-particle orbital 
wave function, and the usual assumption that this single particle orbital wave function is independent of spin. There- 
fore we have i/j afJl (r) = a p (j) a (Y), where a M = a^,b p is the annihilation operator and 4> a is the lowest single-particle 
orbital wave function for species a and spin- independent. Then the Hamiltonian can be simplified as 



c a „n c b 



H=*-S 2 a + -S 2 b + c ab S a ■ S b - 7 B ■ S a - 7 B • S b , (5) 
where a constant is neglected, 



a^F^otv (6) 
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is the total spin operator for species a, c a = d% J d 3 r\4> a \ 4 is the intraspecies spin coupling strength, c ab = 
c% b f d 3 r\(p a (j)b\ 2 is the interspecies spin coupling strength, and we have set -f a = jb = 7, as indeed so for atoms 
with a same nuclear spin. Here we have neglected the quadratic Zeeman effect. This is reasonable under certain 
circumstances, as can be estimated by using parameter values for Na fill ]. The quadratic Zeeman energy is qBg, 
where q — 278Hz/G 2 , Bq can be lOmG to 500mG, hence the quadratic Zeeman energy is about 2.78 x 10~ 2 to 70Hz. 
The linear Zeeman energy is about 1 to 100HZ. Therefore it is easy to reach the regime where the quadratic Zeeman 
effect is negligible. 

S a , Sb together with the total spin S and its z-component S z are all good quantum numbers, as S 2 , S 2 and S 2 all 
commute with the Hamiltonian ([S]). However it should be noted that S a and Sb are not fixed numbers, as in the case 
of pseudospin-i atoms, for which one can find S a = N a /2 and Sb — Nb/2. In the present case, S a , Sb and S should 
all be determined by minimizing the energy. 

In the presence of a magnetic field, for a given S, S z = S minimizes the energy. With S a , Sb, S and S z all being 
good quantum numbers, the ground state is 

\G) = \S™,S?,S m ,S m ), (7) 
where S™, S™ and S m are, respectively, the values of S a , Sb and S that minimize the energy 

„a _ ab b _ „ab 

E = ~^—S a (S a + 1) + — — S b (S b + 1) 

2 r ab 2 (8) 

+—S(S + 1)-^BS, 

under the constraints 

\S a -S b \<S<S a + S b . (9) 

Note that the existence of three quantities S a , Sb and S with the constraint ([§]) as well as the limited ranges of 
S a and Sb, and the dependence on the three parameters c a , c b and c ab ^ makes this minimization problem highly 
nontrivial. We have managed to solve this problem in most of the parameter regimes, as reported in in Appendices. 
Before discussing these cases of c ab ^ 0, we shall first take a look at the case of c ab = 0. 

As we shall discuss different ground states in different regimes of the parameter space, some explanation of the 
nomenclature is in order here. The ground states in two neighboring parameter regimes are said to be continuously 
connected if each of them approaches the ground state on the boundary, when the parameters approach the boundary. 
It is then said that they belong to a same quantum phase. In contrast, if the two ground states in the two neighboring 
regimes approach different limits when the parameters approach the boundary, it is said that there is a discontinuity or 
quantum phase transition. There are several cases of discontinuity, for example, the two limits may be both different 
from the that on boundary, and they may also be two of the degenerate ground states on the boundary, besides, there 
is also the case that the ground state in one of the regimes approaches a ground state on the boundary, while the 
ground state in the other regime approaches a different limit. 

The most interesting ground states in our system are those of EBEC, i.e. BEC with interspecies entanglement. Note 
that throughout this paper, a state which may be entangled is written in the the general form, i.e. S™, S™, S m , S m ). 
A state which is certainly disentangled is written in the form of \S™, S™) a \S™, S™)b 

III. c ab = 

Without spin-exchange interaction between the two species, i.e. c ab = 0, the two species can be considered inde- 
pendently. The ground states are all disentangled. We have E = E a + Eb, with 

E a = ^-Sai^Sa + I) - jBSc , . 

= i(S a -^ + ±) 2 + const, (W) 

for a — a, b. Throughout the paper, we use const to represent a constant whose actual value is not concerned and 
may not be the same each time it appears. 

If c a < 0, we always have S™ — N a . If c Q > 0, one finds the following three subcases, (i) If c Q > 27B, then S™ — 0. 
(ii) If 2N B +\ — ° a — 27-B) then S*™ = n a = Int{^ — |), where Int(x) denotes the integer closest to x while in its 
legitimate range, e.g. < n a < N a . (iii) If c a < 2 n' B +i 1 then S™ — N a . As n a reduces to and N a respectively at 



4 



\N a ,N a ) a \0,0} b 


\n a ,n a } a \0,0) b 


|0,0)a|0 > 0)i 


\N a , N a ) a \n b ,n b ) b 


\n a , n a ) a \n b , n b ) b 


\0,0) a \n b ,n b ) b 


\N a ,N a ) a \N b ,N b ) b 


\n a ,n a } a \N b , N b )b 


\0,0} a \N b ,N b ) b 



FIG. 1: Ground states in c a — c b parameter plane for c ab = and B > 0. They are all direct products of the ground states 
of the two species, and are thus written in the form of S™) a \S™, S b l ) b . The two horizontal lines are c a — 



2N„ 



and 



c a = 2jB, while the two vertical lines are c b = 2 n'^-i an( ^ c& ~ 27B 



Int{ 
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1 ^ 



the two boundaries, the ground state in 



c a < 



2jB 



2-jB 



< c Q < 2^B is continuously connected with those in c a > 2^B and in 



2N a + l ' 

Thus one obtains all the ground states in the parameter subspace of c ab 



b )b and as depicted in FIG.[T] There are nine regimes, each is defined by the range of d 



and B > 0, which can be written in the 

and 



form of \S™,S™) a \S?,S t 

c b specified above. In each regime, the ground state is a direct product of the ground states of the two species given 
above accordingly. Each ground state is continuously connected with those in the neighboring regimes. Therefore, on 
c a — c b plane, ground states in all regimes belong to a same quantum phase. 

As B — > 0, however, the five crossover regimes tend to vanish, and the four ground states in the remaining four 
corner regimes become discontinuous, as already known [32j . Therefore, the quantum phase transitions among the 
ground states in the four quadrants of c a — c b plane for c ab = in absence of a magnetic field can be circumvented by 
turning on and then off a magnetic field. Hence a magnetic field has an interesting effect even in the regime without 
interspecies spin exchange. 

One can imagine the three-dimensional parameter subspace of c ab — 0, with c a , c b and B > as the three 



coordinates. The boundaries ( 
and extending to positive infinities. 



27B 

2N a + l • 



2jB, 



2N b + l 



and c = 27B are all planes starting from the origin 



IV. 



c ab < 



For c ab ^ 0, we have worked out the complicated problem of minimizing E with four variables S a , Sb, S and B in 
most parameter regimes. But in some regimes, the calculations are too difficult or complicated for us to obtain the 
results. The calculation details are given in Appendix A. The ground states we obtained are listed in Table fl] 

For c ab < while B > 0, we have found the ground states in the second, third and fourth quadrants of c a — c b plane, 
as depicted in FIG. [2] In comparison with the case of B = a magnetic field B > both shifts the positions of 
the boundaries and modifies the ground states in the crossover regimes. 

In the three outmost regimes, only the boundaries are shifted, while the ground states remain the same as those 



of B = 0. The ground state is \N a , N a ) a \N b , N b ) b in the regime c a < while c b < 



2-fB-2N a c 
2N b + l 



a < 2~ t B-2N b c a 
— 2N a + l 



In the regimes c a < while c" > 2jB - 2N a c ab and c" < while c a > 2-fB 



- and c 

2N b c ab 



b < while 
the ground 



states are \N a , N a } a \0,0) b and |0, 0) a \N b , N b ) b respectively. 

In the crossover regimes, the ground states are also modified. For c a < while ^| 
the ground state is |JV , iV ) |ni, n'-^b, where n[ = Int 



-2N a c a 



2N b + l 



<c b < 2jB - 2N a ti 



.ab 



-jB+N a \c a 



]. Likewise, for c b < while 



-2N b c a 



2N a + l 



<C a < 



5 



No. 


Parameter regimes 


Ground states 


l,A2a 


c ab < 


c a < 0, 

b ^ 2jB-2N a c ab 
- 2N h + l 


\N a ,N a ) a \N b ,N b ) b , 
disentangled 


A2b 


c a < 0, 

2-yB-2N a c ab ~ ^ „ D r, »t „ab 
2iV,, + l < c b< 27-B 2N a C 


|Wa,A r a)a|ni,ni) () 

disentangled, n\ = Int\ lB+N ^ c — - — il 


A2c 


c a < 0, 

c b > 2yB - 2N a c ab 


\N a ,S z ) a \0,0) b , 
disentangled 


l,A3a 


c" < 0, 

a ^ 2jB-2N b c ab 
- 2JV„+1 


\N a ,N a ) a \N b ,N b ) b , 
disentangled 


A3b 


c" < 0, 

27B 2~N N + h f b < ^ < 2 7 B 2N b C ab 


\n' 2 ,n' 2 ) a \N b ,N b ) b 

j:„„„j-„„„i„j „' — r„j-!l B + N b\ cab \ 11 
disentangled, 712 = inty — 2 J 


A3c 


c" <0, 
c a > 2 7 B - 2jVf,c ab 


\0,0) a \N b ,N b ) b , 
disentangled 


IV 


< c ab < 2-yB 


c a < c ab , c" < c ab 
(boundaries were discussed in Ref. [301) 


\N a ,N b ,n,n),n = Int{^-\) 

(JIl I dueled j 


Bl 


c a >c b > c ab 
N a = N b — N 


c > 27,0 


|0,0> a |0,0> 6) 

disentangled 


B2a 


c 6 < 2jB < $^c» 


\0,0) a \n b ,n b ) b , n b = /nt[^ - |] 
disentangled 


B2b 


27B > $^c b 


1 \ 1 \ r ,\2 1 B\c b -c ab \-c b (c a -c ab )i 

|no3,no3)o|nfti,n 6 i)6, n a3 = Int[ 2[c°- c >>-( c ^) i ] J 

rH"rntin~lrd m.i — fnf H 3-0 * na3 M 


CI 


c a > c aD > c" 

c a c b > ( c abj2 

N a = N b — N > N* 


2 7 B < S^c" 


0,n o ,nt,,nt,) = 0, 0) a n b , n b }(, 
disentangled 


C2 


27B > $¥^c b 


|n a 3, nb2> n b 2 — n a 3, n b2 — n a z) 
entangled, xi b 2 = lTit[— — — b — — — ^] 


Dl 


c b >c a > c ab 
N a = Nb — N 


c a > 2jB 


|0,0>„|0,0>i, 
disentangled 


D2a 


c a < 2jB < $E^c a 


\n a ,n a )a\0, Q) b , n a = Intl^sr - \] 
disentangled 


D2b 


27S > f^c a 


\nai,n a i)a\n b 3,n b 3) b , n b3 - Int[ 1 2 [ e °c''-(c°'') 3 ] J 

diseiiiangieu, iia\ — liny — 


El 


c" > c ab > c a 

c a c b > ( c abj2 

N a = Nb — N > N* 


27B < f^c* 


n a ,0,n a ,n a ) = \n a ,n a ) a \0,0) b 
disentangled 


E2 


27B > £^c a 


\n a 2, n b 3, n a 2 — nbz,n a2 — n b z) 

ill r , t y B4-C ab 7lh^i li 

entangled, n a 2 — Inty- 1 — — — — ~ y\ 


II 


c ab > 2-yB > 


„a ^ „ab Jo ^ „ab 
C > C , C > C 

(boundaries were discussed in Ref. [301) 


|0,0> o |0,0>6, 

disentangled 


III 


c a < c a0 , c" < c a0 , N a = Nb = N 
(boundaries were discussed in Ref. [30]) 


\N,N,0,Q) 
entangled 



TABLE I: Ground states of a mixture of two spin-1 atomic gases in different parameter regimes for c ab / and B > 0. For 
c a < 0, the ground states are always disentangled, and are written in the form of \S™, S?) a \Sr, S b m >b- For c ab > 0, some ground 
states are disentangled and are also written in this form. Some ground states are written in the form of \S™ , S™ , S m , S m ), 
which may be entangled or disentangled, depending on the parameter values. The ordering numbers are just the corresponding 
section numbers in the Appendices where the calculation are done, or as numbered in Ref. [31J for those calculated there (IV, 
II and III). N* is defined in Eq. (^4|) . 



2 7 B - 2N b c ab , the ground state is \n' 2 , n' 2 ) a \N b , N b ) b , where n' 2 = Int[ lB+ ^ c - |]. 

We see continuous connections in both c ab and B dimensions. As B — > 0, all the ground states in these three 
quadrants of c a — c h plane for c ab < reduce to the corresponding ones in absence of a magnetic field. On the other 
hand, as c ab — > 0, the ground states in the these three quadrants reduce to those for c ab = 0, given in last section. 

Note that in all subregimes of c ab < 0, the ground states are always disentangled, as 5"" = S™ + S™. 
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c 

|iVa,iVa>a|0,0) 6 




|A r a,A r a)a|n' 1 ,n' 1 ) 6 


1 A T" TV 7 \ 1 jV T T\ T \ 

\Na,N a ) a \N b ,N b ) b 


\n' 2 ,n' 2 ) a \N b ,N b ) b 


\0,0} a \N b ,N b ) 



FIG. 2: Ground states in c a — c b parameter plane for c ab < and B > 0. They are all disentangled and are thus written 
in the form of \S™ , S^) a \S b n , S^b- The states in the second, third and fourth quadrants are all continuously connected on 



boundaries. n[ = Int[ 



-yB + N a \c a 



C b = 



2~tB-2N a c a 
2JV 6 + 1 



and c b = 2^B 



- 2l> n 2 = 

- 2N a c ab . 



Int 



-,B+N b \c a 



|]. The four boundaries are c a = 27 "^ JV 2J ^ C 



The states in the first quadrant have not yet been determined. 



c a = 2 7 B - 2N b c a 



< c ab < 2-/B 



Now we turn to antiferromagnetic interspecies spin coupling. For < c ah < 2"fB, it has been known previously that 
if c a < c ab andc fc < c ab , the ground state is \N a , N b , n, n), where n = Int{^-\) satisfies \N a -N b \ <n< N a +N b [32lj . 
This state is entangled unless n — N a + N b . The ground states on the boundaries c a 



and c 



have also 



been discussed in details. Especially, it has been known that if c a — c b — c ab , then there are many degenerate ground 
states in the form of \S a , S b , n, n), as far as S a , S b and n satisfy the constraint \S a — S b \ < n < S a + S b . 

We have also determined the ground states in the regime < c ab < 2^B and c a c b > (c ab ) 2 . This regime is 
divided into seven subregimes, but the ground states are all continuously connected on the boundaries between these 
subregimes, as depicted in c a — c b phase diagram for a given value of c ab with < c ab < 2jB (FIG.©, drawn according 
to Table H 

For c b > 2jB and c a > 2jB, the ground state is |0, 0,0,0) = |0, 0) Q |0, 0) b . In the regime defined by c b < 2jB, 



2jB < 



(i.e. below the hyperbola 2jB 



7 c b ), and 2jB < 



7 c b (i.e. above the hyperbola 2<yB 



c c abl c b c b ), the ground state is |0, 0) o |n&, n b ) b , where n b = Int[^3 — \], which reduces to when c b — 2"/B. In the 

regime defined by c a > c b > c ab and 2-fB > ~ ~" ■•" — • 1 ••••« '"•> ••" - 1 u - ' ' 

state is \n a3 , n a3 ) a \n b i, n b \) b , where n a3 = Int{ 



that is, surrounded by c a = c and 2^/B = c b _ c ab c b , the ground 



n b i 



], These states are 



all disentangled. The ground state is \n a3 ,n b 2,n b 2 — n a3 ,n b 2 — n a3 ), where n b 2 = Int[ lB+c cb "° 3 — ^], in the regime 
numbered as C2 and defined by c a > c ab > c b , c a c b > (c ab ) 2 and 2-fB > $^c b , for N a = N b = N > N*, where N* 
is given in (|C4j) . This state is always entangled except on the boundary 2^/B — c c ab~^ c b c b , where the state reduces to 
0, 0} a |n;,, n b ) b . By substituting the boundary values of c a and c b to the values of S a , S b and S that depend on them, 
it is not difficult to see that the ground states in each subregime continuously connected with those in its neighboring 
subregimes. 

Similarly, by exchanging the labels a and b, we know the ground states in the part of c a c b > (c ab ) 2 with c a < c b 
and c a < 2^B. In the subregime defined by c a < 2^B, 2^B < ^„~^ ai , c a (i.e. above the hyperbola 2^/B = ^„~^ at c°), 
and 27B < ^ c a (i.e. on the right of the hyperbola 2~fB = c cab ~^ ca e a ), the ground state is \n a , n a ) a |0, 0} b , 
where n a = Int[^§- — \]. In the regime defined by c b > c a > c ab and 2^B > c a ~ c ab c a , that is, surrounded by 



c a = c b and 2 7 B = 



T c a , the ground state is \n a x, n a x) a \nb3, n b3 )b, where n b3 = Int[ 



2-yB\c a 



2[c<*c l >-(c ab ) 2 ] 
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FIG. 3: Ground states in c a ~ c b parameter plane for < c ab < 2-yB and B > 0, with N a = N b = N . A state which is 
disentangled is written in the form of \S™ , S™) a \S b n , S , b n )b- A state which may be entangled is written in the general form 
of \Sa,Sr,S m ,S m ). The borders of \N,N,n,n) in the lower left part is c a = c ab and c b = c ab . The ground state here is 
\N a , Nf,,n,n) for generic N a and Nb, which is entangled unless n — N a + Nt. The complete hyperbola is c a c = (c ab ) 2 , above 
which all the ground states are continuous connected, forming a single quantum phase. It is |0, 0} a |0, 0)t, in the upper right part 

c a > 2jB and c b > 2jB. Below it the ground state is |0, 0} a \nb, Ub)b, which is bounded on the left by 2yB < c cb Z_ c c ab ^ \ an d on 
the bottom by 2-yB < c ca b~^ c b gt '■ Surrounded by the former boundary and c a = c b is the regime numbered as B2b, where the 

ground state is \n a z, n a z) a \nbi, Tlbi)b- ln the regime numbered as C2, i.e. between 27B < c cab ~^ cb c b and c a c b = (c ab ) 2 , the ground 
state is \n a 3, n b 2, n b 2 — n a3 ,n b 2 — n a3 ), which is entangled except on the boundary with |0, 0) a \n b , n b )b- Similarly, the ground 
state is \n a , n a } a \0, 0} b , in the regime bounded by c a = 2^B on the right, 27B = ^ Q ~° Qb c° below and 2^B = ^1 c a c " on the left. 
In the regime numbered as D2b, as surrounded by c a — c b and 27B = c 'i the ground state is \n a i, n a i) a \nb3, ^3)1,. In 

the regime numbered as E2, i.e. between c a c b = (c ab ) 2 and 2-yB = ^, a b^ c a c " ; the ground state is \n a 2, n^i, n a 2 — ribs, n a 2 — nba), 
which is entangled except on the boundary with \n a , n a ) a |0, 0)&. For the regimes c a > c ab > c b and c b > c ab > c a , the results 
are subject to the condition N > N* . n, n a , rib, n a 2, rib2 n a 3, ribs and A^* are defined in the main text and Table[I] 



n Q i = Int[ lB ^ a " 63 — i]. These states are all disentangled. Finally, the ground state is \n a 2,nb,n a 2 — ni,,n a 2 — n b ), 
where n a 2 = Int[ lB+c b " b —5], in the regime numbered as E2 and defined by c b > c ab > c a , c a c b > (c ab ) 2 and 

b_ ab b_ ab 

27B > c c o.bl c a c a , for N a — N b — N > N*. This state is always entangled except on the boundary 2jB = c c abl c a c a , 
where the state reduces to n ai n a ) a |0, 0)^. The ground states in each subregime continuously connected with those 
in its neighboring subregimes. 

Also, the ground states \n a \, n a i) a \n b 3, n b 3)h and |n a 3, n a 3) a \ribi, n b i)b are continuously connected on the boundary 

c ab K c a = c b < 2 7Bj with = = ^ = = W~£ b) . 

But at the point c a = c b = c ab , the ground states in the six subregimes of c a c b > (c ab ) 2 converging at this point are 
discontinuous with other, as one can see from n a = n b = Int(^ — |), n a 3 = nbi = n b 3 = n a i = Int(^§ b — j) and 
nb2 = n a2 = Int(^^ - |) at this point. 

As c ab — > 2jB, c a = c b = 2jB approaches c a = c = c ab , hence the regimes B2b and D2b tend to vanish. 



VI. c ab > 2<yB 

For c ab > 27B 30], it has been known that the ground state is |0, 0) Q |0, 0) b for c a > c ab and c b > c ab , and 
is |JV, JV, 0,0) for c a < c ab , c b < c ab and N a = N b = N. This result is consistent with with those obtained for 
< c ab < 2-fB. Hence for the regime c a > c ab and c b > c ab and the regime c a < c ab , c b < c ab , N a = Nb — N, the 
ground states for c ab < 2jB and those for c ab > 2jB are continuously connected at c ab — 2jB. 
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|0,0)„|0,0)i 


\N,N,0,0) 





FIG. 4: Ground states in c a — c parameter plane for c ab > 2-fB and B > 0, with N a = Nb = N. The ground state is the 
entangled state \N, N, 0, 0) in the regime c a < c ab and c b < c ab . It is the disentangled state |0, 0} a |0, 0} 6 in the regime c a > c ab 
and c b > c ab . 

Especially, in the cases we have studied, under the condition N a = Nb = N, in the regime c a < c ab and c b < c ab 
and in the regime c a > 2-fB and c b > 2-fB, c ab = 2-fB > is not a boundary, i.e. the ground states are respectively 
the same in these two regimes for c ab > 2-fB and for < c ab < 2-fB. 

VII. QUANTUM PHASE TRANSITIONS 
A. Quantum phase transitions at c a = c b — c ab 

In the regime of < c ab < 2-fB, quantum phase transitions take place at c a = c b = c ab , which is the boundary 
between the two phases discussed above for < c ab < 2-fB. It is a point on c a — c b plane with given c ab and B, and 
is a line in the three dimensional c a — c b — c ab subspace with a given B, and is a two-dimensional surface in the four 
dimensional c a — c b — c ab — B space. 

At c a = c b = c ab , any state in the form of \S a , S b , n, n) with arbitrary legitimate values of S a , Sb and n is a ground 
state. Therefore its degenerate ground state space includes the ground states in all the seven regimes we have studied 
that contact at this degenerate point, that is, the ground states of the six subregimcs of c a c b > (c ab ) 2 neighboring at 
c a = c b = c ab , as well as the ground state \N a , Nb,n,n) in the regime c a < c ab and c b < c ab . Therefore in entering 
c a = c b = c ab from one of the two phases, the ground state remains as the original, and then discontinues in entering 
the any of the other six regimes. Note that there is also a discontinuity in transiting, through the critical point 
c a = c b = c ab , from one of the six regimes belonging to the same phase. 

In any of these seven regimes converging at the point c a = c b = c ab , we always have S m — S™ — n. Therefore, the 
quantum phase transition is a continuous transition. 

If B — > 0, for < c ab < 2"fB, the range of c ab has to be diminished as well, hence in this regime c ab — > too. 
Consequently the regime of |0, 0) a |0, 0)& expands to occupy the whole first quadrant of c a — c b plane, while the other 
regimes in the first quadrant have to be diminished. On the other hand, we have n — > as a consequence of B — > 0. 
Therefore the ground states approach to the corresponding ones for c ab > 2jB, but c ab being infinitesimally positive 
is qualitatively different from the case of c ab = 0, as will be discussed in the next subsection. 

For c ab > 2-fB, the quantum phase transition from the ground state \N, N, 0, 0} in the regime c a < c ab while c b < c ab 
to the ground state |0, 0) a |0, 0)b in the regime c a > c ab and c b > c ab is similar to the one for c ab < 2-fB, with n in the 
latter becoming 0. Hence it is also a continuous quantum phase transition. 
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\n a , n a } a \0, 0)b 


|0,0> a |0,0> 6 




\n a , n a ) a \n b , n b ) b 


|0, 0) a \n b ,n b ) b 


\N a ,N b ,n,n) 





FIG. 5: Ground states in c a — c b parameter plane for c ab = 0+ and B > 0. The two horizontal lines are c a = and c a = 27B, 
while the two vertical lines are c b = and c b = 2jB. The ground states are in the disentangled form of S™) a \S™, S™) b 
in all the regimes with c a > and c b > 0, but are in the form of |S™, S&S 5"™, S m ) in the regime with c a < and c b < 0. 

B. Quantum phase transitions from c ab — to c ab > 

On the other hand, when c ab —> from the positive side, i.e, c ab = 0+, under a given B, the regime of |0, 0) a |0, 0)& 
remains unchanged, while the regime of \N a , Nb, n, n) approaches the third quadrant, and the hyperbola c a c b = (c ab ) 2 
approaches the positive c a and c b axes. The regime of |0, 0} a |nfc, nb)b approaches c a > 2^B while < c b < 2^B, 
similarly the regime of \n a , n a ) a \0, Q)b approaches c b > 2jB while < c a < 2-fB. The regimes B2b plus D2b become 
the regime < c a < 2jB and < c b < 2jB, where the ground state approaches \n a , n a } a \rib, rib}b- Note that the 
subregime of c a c b > (c ab ) 2 with c a > c ab > c b and that with c b > c ab > c a , including the two subregimes of entangled 
ground states, tend to vanish. In FIG. El we draw the c a — c b phase diagram for a given B while c ab — Q+, referring 
to that c ab — > from positive. 

Comparing the ground states of c ab = (FIG. [J) and c ab = 0+ (FIG.O, we can see there are discontinuities between 
c ab > and c ab = 0. First, in the third quadrant, the ground state is \N a , N a )JNb, Nb)b for c ab — 0, discontinuing 
with \N a , Nb,n,n) for c ab > 0. This discontinuity already exists when B = [32|]. This quantum phase transition 
is first order as S m has a discontinuity except in the special case n = N a + Nf,, for which the transition becomes 
continuous. 

Moreover, there are also other discontinuities, which are induced by B > 0. On a c a — c b plane, the boundaries 
c a = 2jB and c b = 2jB exist both for c ab = 0+ and c ab — 0. However, the other two boundaries are different, that 
is, they are c a = and c b = for c ab = 0+, but are c a = 2 ^ B +1 and c b = gffi^ for c ab = 0, though the differences 
diminish as JV a and Nb approach infinities. 

Consequently, there are five discontinuities induced by B for c a > and c b > 0. In the regime < c a < 2 ^ f 7 "^ 1 
and c b > 2-fB, the ground state discontinues from \n a , n a ) a \N b , N b ) for c ab = 0+ to \N a , N a ) \N b , N b ) for c ab = 0. 
This is a first order quantum phase transition except in the special case of n a = N a , in which the transition becomes 
continuous. 

In the regime < c a < 2 ^ f 7 '^ 1 and 9 ^ J ^ 1 < c b < 2jB, the ground state discontinues from |n a , n a ) a \rib, rib) for 
c ab = 0+ to \N a , N a )\nb,rib) for c ab = 0. This is a first order quantum phase transition except in the special case of 
n a = N a . 

In the regime < c a < 2 ^ 7 ^ 1 and < c b < 2N b +i » the ground state discontinues from |n a , n a ) a \rib, rib) for c ab = 0+ 
to \N a , N a ) a \Nb, Nb) for c ab = 0. This is a first order quantum phase transition except in the special case of n a = N a 
while rib — Nb. 

In the regime 2 n <I +i — c ° — an( i < c b < 2 N f ,!+i ' ^ ne g rourL d state discontinues from \n a , n a ) a \rib, nb) for 
c ab = 0+ to \n a , n a ) a \N bl N b ) for c ab — 0. This is a first order quantum phase transition except in the special case of 
n b = N b . 
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In the regime c a > 2jB and < c b < 2N b +i > the ground state discontinues from |0, 0) o |nfe, rib) for c ab = 0+ to 
0, 0) a \Nb, N b ) for c ab = 0. This is a first order quantum phase transition except in the special case of n b = N b . 

Therefore, we find five places of quantum phase transitions from B = to B > 0. In other words, the entire 
subspace of B = is critical. 

VIII. INTERSPECIES ENTANGLEMENT 

Our results indicate that a necessary condition for the ground state to be entangled between the two species is 
c ab > 0. We have found that the ground state is an entangled state entangled \N a , Nb, n, n) for < c ab < 2^B, 
c a < c ab and c b < c ab . In case N a = Nb — N, we have also found that the ground state is a maximal entangled state 
\N, N, 0, 0) for c ab > 2 7 B, c a < c ab and c b < c ab . 

With interspecies entanglement, the occupation number of each spin state of each species is subject to fluctua- 
tion (28|. However, even in absence of interspecies entanglement, such fluctuations can still exist, and there can be 
occupation number entanglement among different single particle states defined by the spin and the species. Such is 
the singlet ground state of single species of spinor atoms, for example. Therefore particle number fluctuations are not 
satisfactory characterizations of interspecies entanglement caused by interspecies spin exchanges. 

A better characterization is an interspecies correlation function, e.g. (N aa N ba i) — (N aa ) (Nba') , which vanishes for 
disentangled state and is nonvanishing if there is interspecies entanglement [28| . 

One can also simply use the spin of freedom of the two species to discuss the entanglement between the two species, 
treating the two species like two giant spins. Then, of course, the entanglement entropy can be calculated. For state 
1 S*™ , S" b m , S m , S m ) , the entanglement entropy is 

£ = - \g(S bz )\ 2 log 2Sr+1 \g(S bz )\ 2 , (11) 

where it is assumed that S™ > Sf, g{S bz ) = (S™,S m - S bz ; S^, S bz \S™ , S£™ , S m , S m ) is the Clebsch-Gordan coef- 
ficient. If < S™, the subscripts a and b are exchanged. £ = for disentangled states, while £ = 1 for state 
\N,N, 0,0). 

We also note that there is a simple yet experimentally measurable quantity as a characterization of the interspecies 
entanglement. This is just the total magnetization S m . If S m = S™ + 5™, there can only be one term in the 
Schmidt decomposition of the ground state in terms of |S'™,5' az } and 1 , Sf, z } , consequently it is disentangled. If 
S rn = \S™ — 5™ |, the ground state is entangled, as there is 2L + 1 terms in the Schmidt decomposition, where L 
represents represents the smaller one of S™ and S 7 b n . 

IX. SUMMARY 

We have obtained most of the ground states of a mixture of spin-1 Bose gases with interspecies spin coupling in 
presence of a magnetic field. For c ab < 0, the ground states, which are all disentangled, belong to a single quantum 
phase. For c ab < 0, a magnetic field modifies the ground states and the boundaries between them. For c ab = 0, a 
magnetic field induces some crossover regimes, hence discontinuities between ground states in different quadrants of 
c a — c b plane in absence of a magnetic field now disappear. 

For c a < c ab and c b < c ab , a magnetic field divides the regime of c ab > into two regimes continuously connecting at 
c ab = 2^/B. For < c ab < 2jB, in this regime of c a and c b , the ground state is \N a , N b , n, n), where n = Int(^ — i) 
satisfies \N a - N b \ < n < N a + N b . For N a = N b , it is continuously connected with \N, N, 0, 0) for c ab > 2~/B in the 
same ranges of c a and c b . It is discontinuous with \N a , N a ) a \N b , N b ) b for c ab — in the same ranges of c a and c b , as 
in the case without a magnetic field. This is a first order quantum phase transition except n = N a + N b . 

Moreover, a magnetic field causes discontinuities of ground states between c ab > and c ab = in the first quadrant 
of c a — c b plane. These discontinuities do not exist in absence of a magnetic field. As c ab — > but remains positive, 
a magnetic field causes the division of the first quadrant into four regimes with continuous connecting ground states, 
as shown in FIG. [SJ while for c ab — there are nine regimes with continuous connecting ground states, as shown 
in FIG. [1] The boundaries of the ground states in these two cases do not match, leading to discontinuities between 
c ab > and c ab — 0. These are usually first order quantum phase transitions except in some special cases. 

c a = c b = c ab > is extremely interesting place, where continuous quantum phase transitions take place no matter 
whether there is a magnetic field and no matter what is the actual value. 
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In terms of bosonic degrees of freedom, the general expression and its composite structure of \S a , Sb, S, S) have 
been discussed previously [HI, H3 ■ ^ wm be very appealing to study the different physical consequences and the 
experimental probes of the crossovers and the discontinuities or quantum phase transitions of the ground states, and 
the effects of interspecies entanglement. 
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Appendix A: ST, ST and S m for c ab < 0, B > 

In this appendix, we find out S™,S™ and S m , in which E is minimal, in the case of c ab < and B > 0. In 
the discussions, E always represent the energy as low as can be determined in the regime under discussion, i.e. the 
meaning of E keeps updating. 

With c ab < 0, E is minimal when S = S a + Sb- Hence the ground state with S z = S is always disentangled. Now 

E = ^S a {S a + 1) + jS b (S b + 1) + c ab S a S b - iB{S a + S b ). (Al) 



Thus 

dE 

dS~ a ~ " " a ' " ~° 2 



c a S a + c ab S b + ^---yB, (A2) 



dE 
dSt 

We consider three subcases in the following 



c b 



( . s , c b S b + c ab S a + - ~ 1 B. (A3) 



1. c a < 0, c b < 



In this subcase, §§- < 0, g < 0, hence S*™ = N a , S*™ = N b , S m = N a + N b . 

2. c a < 0, c b > 

In this subcase, MS- < 0, hence S™ = N a , 



c b 

E(N a , S b ) = jS b (S b + 1) + c ab N a S b - 1 BS b + const. (A4) 



We represent all the values of S a and Sb as points (S a ,Sb) within the rectangular defined by < S a < N a and 
< Sb < Nb on S a -Sb plane (FIG. H]). = defines a stationary line. The points above this line satisfy > 0, 

while the points below the line satisfy j^S- < 0. One can see three possibilities 



< c b < 2 ^ B - 2JV ° e ° 
a. u <. c ^ 2JV,,+i 



The stationary line, depicted as the dashed line in FIG. [51 crosses with the line Sb = Nb- Hence all points with 
S a = N a satisfy |f < 0. Consequently = N a , = N b , S m = N a + N b . Note that this regime so defined can be 
combined with case I A 11 with the same result. 
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FIG. 6: All possible values of (S a , S&) are within the rectangular < S a < N a and < S b < A*';,. = is represented as the 
dashed line in case c a < and < c b < 27 "^^ V 2J ^° C — , and is represented as the solid line in case c a < and 2 ~' B ^ N ° C — < c b < 
2-yB - 2N a c ab . 



2N b + l 



2~fB~2N a c a 
2N b + l 



<c b < 2 7 B - 2N a c ab 



The stationary line, depicted as the solid line in FIG. El crosses with the line S a = N a . The crossing point gives 
the minimal energy. Hence S™ = N a , S™ = m, with 

[ jB + N a \c ab 1 
~ b 2 J ' 



n'i = Int[- 



(A5) 



where Int(x) represents the integer closest to x and in the legitimate range of 5b, i.e. now < Int(x) < AV 
S m = N a + ni . 



c. c b > 2jB - 2N a c a 



All points (S a , S b ) in the rectangular satisfy §§- > Q. Therefore = 0, 5™ = S m = AT 



3. c a > 0, c b < 

One simply exchanges the subscripts or superscripts a and b in the preceding subcase. Thus there are also three 
possibilities. 



a. < c a < 



2-,B-2N b c a 
2N a +l 



S™ = N a , S™ = Nt, S m — N a + Nb. This regime can be combined with case I A 11 without the same result. 



2-,B~2N b c 



2N„ 



j— <c a < 2-/B - 2N b c a 



S™=n b2 , with 



n' 2 = Int[ 



jB + N b \c ab \ 1 



(A6) 



now < Int(x) < N a , S™ = N b , S m = n 2 + N b . 



c. c a > 2-yB - 2N b c a 



S™ = 0, S r b n = S m = N b . 
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FIG. 7: For < c ab < 2jB and N a = N b = N, the whole region of (S a , S b ), which satisfy < S a < N and < S b < N, can be 
divided into four regions. (I) St — S a > So > 0, (III) \Sb — S a \ < So < Sb + S a , and the rest regions of (II) and (III). The dashed 
line represents = in the case c a > c b > c ab and c b < 2jB. The solid line represents J^P- = in the case c a > c ab > c b and 
c a c b > ( c abj2_ 



Appendix B: 5™, 5™ and 5"" for < c ab < 2 7 B, c a > c 6 > c a6 , N a — Nb = N 

Define 



So = ^-i (Bl) 



which is the value of 5 on which S'-dependent part of E is minimal if there were no constraint on S. 

With < c ab < 2 7 B, it can be found that the whole region of (S a , Sb) can be divided into four regions, as shown 
in FIG. El 

In region I, Sb — S a > So > 0, hence E is minimal when S = Sb — S a , with 



E = ^-S a (S a + 1) + jS h (S b + 1) - c ab S a S b - c ab S a - 1 B(S b - S a ), (B2) 

for which it is found that J^- > 0. Thus in region I, E reaches its minimum at Sb = S a + So. It is then easy to note 
that the minimum of E in this region rests on S a = 0, Sb = So. 

Similarly, it can be shown that in region II the minimum of E rests on Sb = 0, S a = So- Since both (0,6*0) an d 
(So, 0) also belong to region III, the minimum of E in the whole rectangular must be in regions III and IV. 

In region III, \Sb — S a \ < So < Sb + S a , hence S can reach So, hence 

E = Ca ~ C ° S a (S a + 1) + ^—Sb(S b + 1) + ^yS (S + !) - 7SSo, (B3) 

for which j^jr- > 0, ^j- = c ~ c (2^ + 1) > 0. Thus the minimum of E in region III must rest on the border between 
III and IV, defined by S a + S b = S . 

To conclude the above discussion, the minimum E must locate in region IV, where S a + Sb < So, hence the minimum 
of E lies on S = S a + S b , 

E = ^S a [S a + 1) + jSb(S b + 1) + c ab S a S b - iB{S a + S b ). (B4) 

One obtains 

OE c a 

— = c a S a + c ab S b + — - jB, (B5) 
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dE 
dS~ b 



c Sb + c a S a 



jB, 



(B6) 



according to which one needs to consider two subcases. 



1. c b > 2-/B 

In this parameter regime, W- > 0, 4§- > 0. Therefore E is minimal when S™ = S™ = S m = 0. 



2. c b < 2 7 B 



In this parameter regime, 4§- = defines a stationary line, shown as the dashed line in FIG. [JJ Consequently, the 



minima of E in different parts of region IV are 



E = { 2cb "a 1 l 2 

1 £S a (S a + l)--fBS a , 



'-jB]S a 



if < S a < 2# 



if 



2r" 



■ 2c ' 

< S„ < S , 



(B7) 



For S" a > 2§ — 2§3F! ^ is found that > 0, hence i? reaches its minimum at — 0), i.e. the point bordering 
the other part of region IV. 

Therefore, the minimum of E in the whole rectangular must locate on the dashed line in sector IV, on which E is 
given by (|B7I) . Then there are two possibilities. 



c b < 2yB < 



We have ST 



n qm qm 'jB 1 

' b — — 2- 



b. 2-yB > 



One finds that 



ST = Int 



-fB(c b - c ab ) - c b (c a - c ab )/2) 

c a c b _ ( c afc)2 



°6 — 



1 B - c ab S a 1 
c b 2' 



(B8) 

(B9) 

(B10) 
(Bll) 



Appendix C: S^ 1 , 5^ and S m for < c a6 < 2 7 B, c a > c ab > c b , c a c b > (c ab f, N a = N b = N 



Again, we use FIG. [71 It can be shown that in region IV, < 0, thus the minimum of E in this region lies on the 
border line with region III, i.e. S a + Sb = Sq. It can be shown that the minimum of E in region II lies on the border 
line S a — Sb = Sq with region III. Therefore we need only to consider regions I and III. 

In region III, as shown in last section, S m — Sq, and E is given in Eq. (|B4|) . But now that c < c ab < c a , we have 
J^P- > and j^r- < 0. Consequently the minimum of E lies in the border line Sb — S a = Sq. 

Therefore, we conclude that E takes its global minimum in region I, where, as discussed in last section, S m = Sb — S a , 
E is as given in Eq. (|B2j) . for which 



dE 
dS~ b 



b c ab o 
C bb~ C S a 



■yB. 



As shown in FIG. §^ 
of E are found to be: 



defines a stationary line which crosses with Sb — N at ( 



c b N+c b /2—yB 



(CI) 

, N) . The minima 



15 



E = { ^f^Sl + - ^B]S a + + - I), if < S a < ^±0^, 

In the second interval c JV+ ^ { 2 ~ 7B < S a < N — So, 

— = c a S a ~c ab (N+l) +1 B + - 



a b _ f ab\2 a b 

> r J -N+—(—r + l) 

ab 2 X C ah 



7^ ~ 1). (03) 



which is positive if N > N* , where 



jB(c a - c ab ) + <c ab ) 2 - c ab c a - 
N* = — ' , V / , N9 2 —. (C4) 

Then the minimum of E must locate on the stationary line c b Sb — c ab S a + ^ — 7 £> = 0, with < S a < c N+c r J^~ lB ■ 
One can see two possibilities. 

1. 2 7 B<^<c 6 



In this case, ^f- - ^—^jB > 0. Thus we have 5™ = 0, ST = S m = Int[2% - §] 



2. 2 7 B > 



,6 



Then 



fl7f r 7j?( C ab -c b )-c"(e a -c^)/2) 



'vB + c ab S m 1 

ST = M 7 J a -5], (06) 

gm = s?-S™. (C7) 



Appendix D: S™, ST and S m for < c ab < 2 7 B, c b > c a > c ab , N a = N b = N 
By exchanging the labels a and b in Appendix [B] one obtains the following results. 

1. c a > 2yB 

ST = ST = s m = o. 

2. c a < 2jB 

ST' 



We have ST = 0, 5™ - 5 m = ^ - §. 
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an 

S" 



Int 



-fB(c a - c ab ) - c a (c b ~ c ab )/2) 

c a c b _ ( c afc)2 



-fB - C ab S b 1 

~ a 2' 
ST + SU 1 . 



(Dl) 

(D2) 

(D3) 
(D4) 



Appendix E: S?, 3? and S m for < c ab < 2-/B, c b > c ab > c a , c a c b > {c ab f , N a = N b 
By exchanging the labels a and b in Appendix [C] one obtains the following results. 



= 0, S™ = S m = Int[$r - i]. 



571 

b 



Int[ 



Int[ 



jB{c ab - c a ) - c a (c b - c ab )/2) 

c a c b _ ( c afc)2 

~/B + c ab ST 1 



= sr 



(El) 

(E2) 
(E3) 
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